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ABSTRACT 


The AC impedance spectra of thin (0.27 fim) polypyrrole films have been 
obtained at open circuit potentials from -0.4 V to 0.4 V vs. SCE. Two 
limiting cases are discussed for which simplified equivalent circuits are 
applicable. At very positive potentials, the predominently nonfaradaic AC 
impedance of polypyrrole is very similar to that observed previously for 
finite porous metallic films. Modeling of these data with the appropriate 
equivalent circuit permits effective pore diameter and pore number densities 
of the oxidized film to be estimated. 

At potentials from -0.4 V to -0.3 V, the polypyrrole film is 
essentially nonelectronically conductive and diffusion of polymer oxidized 
sites with their associated counterions can be assumed to be linear from the 
film/substrate electrode interface. The equivalent circuit for the 
polypyrrole film at these potentials is that previously described for metal 
oxide, lithium intercalation thin films. Using this model, counterion 
diffusion coefficents are determined for both semi-infinite and finite 
diffusion domains. The diffusion coefficient data so obtained is consistent 
with the existence of a gradient in the morphology of the pol 3 rpyrrole films 
in a direction normal to the surface of the electrode . 

In addition, the limiting low frequency resistance and capacitance of 
the polypyrrole thin films is determined and compared to that obtained 
previously for thicker films of the polymer. The origin of the observed 
potential dependence of these low frequency circuit components is discussed. 



INTRODUCTION 


A number of technological applications for electronically conductive 
polyheterocycles (eg. polypyrrole, polythiophene) have been demonstrated in 
recent years. These include the application of polyheterocyclic films as 
potentiostatically controlled ion-gates (1,2), in electrochromic displays 
(of. 3, 4), as battery cathode materials (5-8), and in chemically responsive 
transistors (9-11). All of these applications rely on the ability of the 
electronically conductive polymer to "switch" rapidly and reversibly between 
insulating and conducting redox states. The rate of the polymer redox 
reaction for electrochemically synthesized conducting polymers has been 
shown to be inversely related to the thickness of the polymer film (8,10,12) 
For this reason, the design of pol 3 nner-based devices has emphasized the use 
of extremely thin (1 < 1 ^m) polymer films (8,10,12). 

Despite the importance of polyheterocyclic films of submicron 
thickness, little is currently known about mass transport in these thin 
films. Previous AC impedance studies of pol 3 rpyrrole by Burgemayer and 
Murray (1,2) and more recently by Tanguy et al. (13,14) have employed 
relatively thick (> 5 nm) pol 3 nner films. Such films cannot be 
quantitatively addressed electrochemically; ie. , as -synthesized oxidized 
polymer cannot be quantitatively reduced (15,16). Consequently, the 
extraetion of mass transport and kinetic information from this data is not 
straightforward. Bull et al. examined the AC impedance behavior of thin 
polypyrrole films over a broad frequency band but did not consider mass 
transport in the analysis of these data (17). Glarum et al. has recently 
published an analysis of the potential dependant AC impedance behavior of 
thin electronically conductive polyaniline films (18). 
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Recently, AC impedance analyses have been applied to the determination 
of diffusion coefficients in electrochroraic WO 3 films (19) and in redox 
polymer film modified electrodes (20). This work has demonstrated that 
transport information can be obtained from AC impedance data for any system 
which exhibits finite, linear diffusion, providing the heterogeneous 
electron transfer kinetics of the system are suitably fast. AC impedance 
possesses several advantages for the determination of diffusion coefficients 
as compared to conventional, large amplitude, electrochemical experiments 
such as chronocoulometry. In contrast to Cottrell- type experiments in 
which diffusion coefficient information is obtained while the system is in 
flux from one redox state to another, the perturbation of the film redox 
state during an AC impedance experiment is extremely small permitting 
transport information to be obtained for any discrete ratio [Ox]/[Red] . 
Consequently, morphological changes in the film structure which can 
accompany a bulk conversion from one redox state to another are avoided; 
characteristics of the film associated with a particular redox state can be 
elucidated. In addition, a wide band of frequencies can be addressed in an 
AC impedance experiment allowing the characteristic time -scale of an 
electrochemical process, such as diffusion, to be isolated as the 
corresponding frequency interval and differentiated from 

characteristic time -scales. Consequently, both kinetic and mass transport 
information can obtained from the same experiment. 

In this paper, we interpret the AC impedance data obtained for 0.27 nm 
± 0.02 fim films in the context of the limiting behavior which is observed 
for polypyrrole films at low doping levels and for fully oxidized films. In 
the former case (Limiting Case I) , the polymer film is essentially 
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nonelectronically conductive. The equivalent circuit appropriate for this 
system is that described by of Ho et al. (19) and by Franceschetti et al . 
(21-23). As noted above, diffusion coefficient information has been 
obtained using this technique for a number of functionally similar systems 
including WO3 lithium intercalation films (19,21-23) and Nafion polymer 
modified electrodes (20,24). 

The AC impedance of oxidized polypyrrole (Limiting Case II) has 
previously been described in terms of a porous electrode model by Bull et 
al. (16) and by Burgmayer et al. (1,2). However, a more rigorous model has 
recently been demonstrated for finite, porous electrodes by Reiser et al. 
(25) and by Candy et al. (26). The latter model permits film parameters 
such as the equivalent cylindrical pore diameter and pore number density to 
be calculated from the AC Impedance spectra of finite porous films (26) . In 
the present case , the AC impedance observed for oxidized polypyrrole films 
is interpretted using the analysis of Candy and coworkers (26) . 

In addition, the limiting low frequency resistance and capacitance of 
the polypyrrole thin films is determined and compared to that obtained 
previously for thicker films of the polymer (13,14). The origin of the 
observed potential dependence of these low frequency circuit components is 
discussed. 

THEORY 

Limiting Case I: Lightly Doped Polymer. 

As noted above, a number of recent papers describe the interpretation 
of AC impedance data for transition metal intercalation compounds (19, 21- 
23) and redox polymer modified electrodes (20,24). In both of these 
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systems, a redox couple is confined to a thin film, and the diffusion of the 
electroactive species in each is finite and linear. We have chosen to use 
the analysis of Ho et al. because it explicitly treats the limiting cases 
applicable to polypyrrole (19) . The authors have made a number of 
assumptions in the derivation of this model which need to examined here in 
the context of the application of this theory to electronically conductive 
polymers. Two of these assumptions are particularly important: 1) 
diffusion of the intercalating species is linear, and 2) the Warburg 
impedance is that due to the diffusion of a single species (Li°) and hence 
can be characterized by a single diffusion coefficient (19). 

In the case of LiyW 03 and other electronically conductive transition 
metal oxide films , it is commonly assumed that intercalation of Li'*’ is 
preceded by charge transfer and ion pair formation, Liy''TJ 03 ' (or LiyW 03 ) , at 
the oxide film/electrolyte interface (19,22), Because intercalation of the 
Li"*" is commensurate with charge transfer, Li° is the only mobile species in 
the film. Diffusion of the intercalated Li° then occurs linearly from the 
film/electrolyte interface. A key point regarding the applicability of the 
data analysis technique employed by Ho et al. to polyheterocyclic films is 
the extent to which diffusion in the conducting polymer film is linear. At 
potentials where polypyrrole is a good electronic conductor (Eq^ > 0.3 V vs. 
SCE) , diffusion cannot be assiamed to be linear since the polymer redox 
reaction can proceed from any conductive surface in the porous film. 

However, at open circuit potentials from -0.4 V to -0.3 V (100-200 mV 
negative of E°'), where the polypyrrole is lightly doped, the data from the 
present paper as well as previous experimental evidence (27,28) indicates 
that the conductivity of the polymer is very low. As the bulk electronic 
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conductivity of the pol)mier is low, the redox reaction of the polymer can be 
assumed to proceed only from the Pt substrate/polymer interface. Hence, it 
is assumed here that diffusion of oxidized polymer sites (and charge 
compensating anions) is linear from the Pt substrate/polymer interface at 
potentials < 0.3 V vs. SCE. 

The diffusion coefficient for the diffusion of intercalated Li° in a 
metal oxide matrix is many orders of magnitude smaller than that for the 
diffusion of Li'*’ in the contacting electrolyte (ca. 10*^2 vs. 10"^ cm^ sec'^ 
(19)). In contrast to AC impedance analyses of redox couples in solution 
(29) , the diffusion limited current for LiyW 03 is always dictated by slow 
diffusion of Li® and the expression for the Warburg impedance incorporates a 
single diffusion coefficient (19,23). It is less clear in the case of 
polypyrrole if the Warburg impedance can be assumed to be dominated by a 
single diffusing species, eg. diffusion of the anion. The polymer redox 
reaction for pol 3 rpyrrole can be written: 


PPy° + A' 


-> PPy'*"A' + e‘ 


Diffusion of PPy"*" sites generated by an AC varying potential signal might be 

1 tHlitP-Cl bv t"be 1 n*F artinn A* /I •? -F-Fi i<? ■? r*+ 

electron diffusion; (i.e., the rate of electron hopping), or some 
combination of these. Buttry and coworkers have performed quartz crystal 
microbalance (QCM) gravimetry on thin pol 3 rpyrrole and polythiophene films 
synthesized in acetonitrile (30). These data indicate that bulky 
tetrabutyl ammonium cations are excluded from polypyrrole films during 
oxidation- reduction cycles (30). Under these circumstances, if electron 
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hopping can be assumed to be facile, diffusion of PPy"*" will be limited by 
codiffusion of charge compensating BF 4 ‘ anions. Again, a single diffusion 
coefficient, Dgp^., is adequate to describe the diffusion limited current. 
The application of the data analysis technique of Ho et al. to 
nonpermselective polypyrrole films will yield diffusion coefficients which 
are composite values containing contributions from both anion and cation 
diffusion (20,29). 

AC impedance analysis of polypyrrole, then, can be made to conform to 
the assumptions made by Ho et al. by insuring that the polymer is not 
electronically conductive at the potentials at which diffusion coefficient 
information is obtained, and by employing a tetrabutylammonium cation which 
is excluded from the polymer film. 

The model proposed by Ho et al. assumes the Randles equivalent circuit 
shown in Figure lA (31). However, the usual form of the Warburg 
(diffusional) impedance is modified to compensate for the finite diffusion 
length of ions in a thin film. The detailed derivation of these expressions 
is found in reference 19. The explicit form of the finite Warburg impedance 
is (19): 


where : 


- Rw 


j^w 


Xw - (wC„)-l - |Zw|sin /3 
Rw - |Z„|cos /9 
|Zw| - |i^o/io^l 


( 1 ) 

( 2 ) 

(3) 

(4) 


where i/^ is the AC potential amplitude, a is the electrode area, and C^,, and 
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Rw are the Warburg capacitance and resistance respectively. The Warburg 
phase angle, / 9 , is defined as: 


(sinh(2kl))+(sin(2kl)) 

- arc tan (5) 

(sinh(2kl)).(sin(2kl)) 


with: 

k - (w/2D)^/2 

where <ji> is the angular frequency of the applied AC signal and D is the 
diffusion coefficient of the counterion, BF4*. The current amplitude, 
is defined as : 


zF dCgF*. h 2 + s 2 1/2 

io - Vq (wD) 1/2 (y) 

2 dE d 2 


where dCgy^./dE is the gradient of polymer oxidation state, and hence BF4" 
concentration, with open circuit potential. The expressions for the 
equivalent series resistance, Rg, and capacitance, Cg, elements of the 
faradaic impedance (incorporating the finite Warburg elements (Eq. 2 , 3 )) are 
given by ( 29 ) : 


R-S “ ^■ct ®W “ ^ct I *^0/^0® I P 
Cs ” Cw - sin 


( 8 ) 

( 9 ) 
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where is the charge transfer resistance; R^t “ RT/zFIq. The expressions 
for the real, Z[re] , and imaginary, Z[im], components of the total impedance 
are obtained by substituting Eq. 8 and Eq. 9 into the expressions for Z[re] 
and Z[im] in the Randles circuit (29,30): 


Rs 

Z[re] = Ru + - (10) 

((Cd/Cs)+1)2 + (wRgCd)^ 


wRs2c^ + Cd/(Cs2o,l) + («Cs)-l 

zti™] 3 3 (11) 

((Cd/Cs)+1)2 + (wRsCd)2 


Where is the series resistance of the solution and the polymer film. 

With these expressions, Nyquist diagrams can be generated for any 
combination of experimentally adjustable parameters. In practice, two 
important limiting subcases of Limiting Case I can be differentiated: lA) 
thick films and/or small diffusion coefficients (kl » 1) , and, IB) thin 
films and/or large diffusion coefficients (kl « 1) (19). 

Limiting Case lA corresponds to semi-infinite linear diffusion; ie., 
the maximum excursion of the diffusion layer at steady- state is less than 
the film thickness. In this case, becomes 45® and the Nyquist plot (Z[im] 
vs. Z[re]) is linear with a slope of 1. From this linear region, if dE/dC 
is known, the diffusion coefficient can be evaluated from the linear 
dependence of the impedance modulus, |Z| , with (19): 

|Z| - 1 [(dE)/(dC)] [FDi/2a]-Vi/2 I (12) 
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In Limiting Case IB, the diffusion layer accesses the entire, film 
volume corresponding to thin- layer electrochemical behavior. The Warburg 
phase angle, is 90® and the Nyquist plot becomes vertical at a resistance 
of (R]^+R^^+R^) . The value of R]^, and hence the position of this vertical 
branch, is sensitive to the diffusion coefficient (19): 

1 

Z[re] - Ri - (dE/dC) (13) 

3nFaD 

The capacitance also becomes frequency independent under these conditions 
and equal to (19): 

1 

Z[im] - - |(dE/dC)/nFwla| (14) 

cjC]^ 

Eqs . 13 and 14 may be used to calculate D and dC/dE simultaneously from R]^ 
and C]^ obtained from the low frequency region of the Nyquist plot. 

For perfectly homogeneous films, the diffusion coefficients calculated 
from the semi -infinite frequency regime (Limiting Case lA) ought to agree 
with those obtained at frequencies corresponding to the finite diffusion 
case (Limiting Case IB) . 

Limiting Case II: Oxidized Polymer. 

Although the theory for the AC impedance of porous electrodes is well 
established (32-34), only relatively recently has the AC impedance of thin, 
porous, metallic films been examined (25,26). The work of Kaiser et al. 

(25) and Candy et al. (26) has involved finite porous electrodes in contact 
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with electrolyte solutions containing no redox couple; systems for which the 
impedance is entirely nonf aradaic . Candy et al. have shown that porous 
metal layers in which the pores have irregular shapes (as, for example, in 
Raney-metal surfaces), exhibit impedance behavior similar to that predicted 
for surfaces consisting of cylindrically shaped pores of a uniform diameter 
(26) . Such poorly defined metal films can be evaluated in terms of the pore 
radius, pore number density, and pore length of the equivalent cylindrical 
pore electrode (26) . 

Feldburg (35) and Burgemayer et al. (2) have calculated that the bulk 
capacitances of 100-200 C cm‘^ observed for oxidized polypyrrole films can 
only be accounted for if the specific capacitance of the polymer, and hence 
its conductivity, is very similar to that of metals. Oxidized polypyrrole 

films, then, might be expected to exhibit AC impedance behavior similar to 
that of porous metal layers. We have applied the recent model proposed by 
Candy et al. (26) to the analysis of the AC impedance observed for oxidized 
(^oc > +0.1 V vs. SCE) polypyrrole films in electrolyte solutions containing 
no electroactive species. The assumptions inherent in this analysis are 
simply that the conductivity of electrolyte in the pores of the electrode is 
known and is the same in all pores , and that the specific capacitance of the 
porous electrode surface, is known and uniform. In addition, the 

capacitance of the Pt substrate electrode surface is assumed to be 
negligibly small relative to that of the porous bulk of the polymer. 

The equivalent circuit for a single, electrolyte-filled cylindrical 
pore is shown in Figure IB. The cylindrical pore electrode described above 
contains n such pores connected in parallel (26) . This is the equivalent 
circuit describing an equal number of finite, parallel transmission lines 
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(25,34). Note that the faradaic branch found in the Randles circuit 
(containing R^t and Z„) is absent from this circuit. The real and imaginary 
components of the impedance for this circuit are given by (25) : 

sinh (A/1) - sin (A/1) 

Z[re] - Ru + (A/1) - - (15) 

cosh (A/1) - cos (A/1) 


sinh (A/1) + sin (A/1) 

Z[im] - (A/1) (16) 

cosh (A/1) - cos (A/1) 

where 1 is the equivalent cylindrical pore length, and A is the penetration 
depth of the AC signal, given by (26): 

A - l/2(rc<r/<*>Cdi)^^^ (17) 

Where a is the conductivity of electrolyte in the pores , w is the angular 
frequency, r^, is the equivalent cylindrical pore diameter, and C(ji is the 
specific double -layer capacitance. With these equations, Nyquist diagrams 
can be generated for any desired set of electrode parameters . 

Qualitatively, the Nyquist diagrams calculated from this equation are 
identical in shape to those generated by Eqs. 10 and 11 except for the 
absence of the charge transfer limited semi-circle. That is, a 45° linear 
region is observed at high frequencies (A « 1), and both the capacitance 
and resistance reach the limiting values R^ and at low frequencies (A » 
1) . However, the origin of the limiting behavior is different in the 
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present case than in Case I. The penetration depth of the AC signal (Eq. 

17) into the pores of the metallic film is a function of a, C^i, r^,, and w. 
As the frequency decreases, 1 increases and a greater total pore volume with 
its associated resistance and capacitance are accessed. Thus, as predicted 
by Eqs . 15 and 16, increases in Z[re] and (ti)Z[im])'^ are commensurate with 
decreasing w. At some sufficiently low frequency, A - 1 and the total 
volume accessed by the AC signal approximately equals the total available 
pore volume. At this frequency, Z[re] and (wZ[im])l- assume the limiting 
values and C^. 

If a, and 1 are known, the equivalent cylindrical parameters r^, 

and n of the porous film can be calculated from and C^ (26): 


1 

Rl ; (18) 

3n rc^ 

Cl - 2n rclC(ji (19) 


The pore length, 1, and the volume of electrolyte in the pores of the film, 
Vp, are related by the expression: 


1 - (3RiVpa)V2 


( 20 ) 


As a consequence of this equation, either 1 or Vp is sufficient when 
combined with Cjji, o and the experimental values R1 and Cl to calculate r^ 
and n for a porous film (26). 
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EXPERIMENTAL 


Haterials and Equipment 

Platinum foil (t = 0.254 mm, Alfa) and extruded Kel-F rod (d =■ 0.5 in, 
Afton Plastics Molding) were used to construct Pt disk electrodes as 
follows. The Kel-F electrode body was first machined to the dimensions 
shown in Fig. 2. Pt disks (d — 3.3 mm) were then heat-sealed onto the 
surface of the Kel-F body by melting the Kel-F around the Pt disk with a 
heat gun and physically pressing the Pt disk into the heat-softened surface 
with a glass rod. Electrical contact was made to the back of the Pt disk 
with a copper wire and silver epoxy (Epo-Tek 410E, Epoxy Technology) as 
shown in Fig. 2. Excess Kel-F was removed from the surface of the Pt/Kel-F 
electrode by polishing with 400 grit and 600 grit Carbimet (Beuhler) paper 
until a planar electrode surface was established. This surface was then 
further polished sequentially with 1.0 pm, 0.3 pm, and 0.05 pm alumina 
powder (all from Buehler Inc.). After this initial polishing sequence, the 
electrode surface was polished between experiments with 0.05 jum alumina 
only. 

A large area, Pt gauze counter electrode (25 x 25 mm, AESAR) and a 
conventional saturated calomel reference electrode (SCE) were employed for 
all electrochemical experiments . 

Tetrabutylammonivim tetrafluoroborate electrolyte was prepared from 
ammonium tetrafluoroborate (97-H%, Aldrich) and tetrabutylammoniun bromide 
(99% Aldrich) and recrystalized from ethyl acetate -pentane. The 
recrystallized salt was dried envacuo at 100 C for ca. 24 hrs. prior to use. 
Pyrrole (98%, Aldrich) was distilled in an inert atmosphere and stored over 
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argon. Acetonitrile (Burdick arid Jackson, UV Grade) was used as received. 
All solutions employed for electrochemical measurements were purged with 
purified N 2 prior to use. 

A EG&G PAR Model 273 potentiostat/galvannostat was employed for film 
deposition and all electrochemical measurements. Cyclic voltammograms were 
recorded on a Houston Instruments Model 2000 X-Y recorder. AC impedance 
measurements were performed in conjunction with a computer controled EG&G 
Model 5301 Lockin Amplifier and EG&G AC impedance software. At frequencies 
greater than 5 Hz., the applied sine-wave excitation signal was generated 
and the current response analyzed by the lockin amplifier-controlled PAR 
273. The measured impedances were then transfered to the computer for 
display and storage. At lower frequencies, the EG&G software removes the 
lockin amplifier from the circuit and the sine wave excitation signal is 
computer generated. The current response obtained from any desired number 
of cycles is measured by the 273, then analyzed by computer via a fast 
fourier transform (FFT) algorithm. At these lower frequencies, the 
impedance was determined from the data obtained from 3 cycles. The 
amplitude of the applied potential sine wave in each case was 10 mV. 
Calibration of the 5301 is accomplished automatically with internal 
standards. All impedance measurements were normalized to an electrode area 
of 1 cm^ . 

The conductivity of 0.2 M BU 4 NBF 4 , MeCN electrolyte was measured with a 
Yellow Springs Instruments Model 31 AC conductivity bridge and a 
conventional YSI parallel plate conductivity cell. The cell constant was 
determined with KCL solutions . All solutions were equilibrated in a 25 C 
water bath. 
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Film Deposition and Electrochemical Measurments 

Polypyrrole films were prepared from monomer solutions containing 0.5 M 
pyrrole, 0.2 M Bu^NBF^ in acetonitrile. Polymerization was accomplished 
galvanostatically in a conventional one -compartment cell at a current 
density of 1.0 mA cm‘^. A polymerization charge of 0.10 C cm*^ was found to 
correspond to an oxidized film thickness of 0.27 pm ± 0.02 pm as measured by 
a Tencor Alfa- Step profilometer . After deposition, the integrity of each 
polymer film was ascertained by recording a cyclic voltammogram at 20 
mV/sec. All electrochemical measurments were conducted in N 2 purged, 0.2 M 
BU 4 NBF 4 , acetonitrile electrolyte. 

AC impedance data was obtained by first potentiostating the polymer at 
the desired potential for 60 seconds, then equilibrating the polymer at open 
circuit for 120 seconds. The polymer was then potentiostated at the 
terminal open circuit potential before AC impedance data acquisition was 
initiated. This procedure resulted in Nyquist diagrams which exhibited no 
hysteresis with varying potential. This indicates that the impedance data 
is obtained at each potential on polymer which is in equilibrium with the 
mean applied potential. The AC impedance was measured for each open circuit 
potential at 42 frequencies from 0.05 Hz. to 100 KHz. 

RESULTS AND DISCUSSION 

Figure 3 shows a cyclic voltammogram in 0.2 M BU 4 NBF 4 , MeCN electrolyte 
for a typical 0.27 pm pol 3 rpyrrole film employed for the AC impedance studies 
below. The large capacitive currents in the region positive of the polymer 
redox reaction (E > +0.1 V) indicate that the polymer is electronically 
conductive in its oxidized state. At potentials negative of ca. -0.3 V, the 
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capacitive currents decrease dramatically since the reduced polymer is an 
electronic insulator. The unusual, asymmetric shape of this voltammogram 
has been observed previously for thin polypyrrole films prepared and cycled 
in BU 4 NBF 4 , MeCN electrolyte (36). As mentioned above, such cyclic 
voltammograms were routinely used to ascertain the quality of the 
polypyrrole films sythesized prior to performing the AC impedance 
experiments . 

Oxidized polypyrrole 

Theoretical Nyquist data for porous metallic films can be calculated 
for the film parameters r^ and n from Eqs. 15 and 16 provided that a, C^i, 
and either 1 or Vp are known. In the case of porous metal films, values for 
Vp are conveniently obtained gravimetrically from measurements of the water 
content of water impregnated films (26). This procedure, however, is 
impractical for the 0.27 pm polypyrrole films employed in this study since 
such films are too thin to be removed intact from the substrate electrode 
surface. For the purposes of estimating the parameters rj, and n, then, Vp 
was taken to be equal to the total film volume, Vf - r^l, and 1 was 
calculated from Eq. 20. An estimated value for of 20 pF cm^ was also 
employed based on the calculations of Feldburg (35). Again, this number 
which is typical of metal/electrolyte interfaces should be considered to be 
a limitingly high value. Because r^ is proportional to both and Vp, 
values for r^ calculated from these estimates of and Vp should be 
considered maximum values. 

Experimental and simulated Nyquist diagrams for polypyrrole at three 
open circuit potentials in the interval from +0.4 to +0.2 V are shown in 
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Figure 4. Figure 4A shows a comparison of the experimental data obtained at 
low frequencies (f < 5 Hz.)- Note that at all three potentials, Z[re] is 
nearly frequency independent. As shown in Figure 5, plots of Z[im] (- 
l/<*>Cs) vs. for this low frequency data are linear (r > 0.999) 

indicating that the capacitance is constant (Cg - C^) . Because this 
potential region is well positive of the polymer electrochemistry (E° - -0.2 
V) , changes in the formal oxidation state of the polymer over this interval 
are small, ie. the polymer is essentially 100% oxidized at Eqc > 0.2 V. The 
similarity of the three plots shown in Figure 4A, then, corroborates the 
conclusion that the observed impedence is predominately nonfaradiac. Thus, 
Case II limiting behavior is observed. 

Figures 4B - 4D show the high frequency Nyquist data in greater detail. 
Simulated plots were generated in each case by fitting to the experimental 
data using r^, and n as adjustable parameters. The resulting best-fit 
values for these parameters are listed in Table I. Note that the agreement 
between the simulated and the experimental data is best at 0.393 V (Fig. 

4B) . In particular, at this potential the elbow in the Nyquist plot occurs 
at very similar frequencies; ca. 40 rad sec‘^ (expt.) vs. 25 rad sec"^ 
(calc.). The values for rj, and n of 4.26 A and 5.43 x 10^^ pores cm‘^ 
obtained from the simulation at 0.393 V (Table I) are reasonable limiting 
values in view of what is known about the morphology and mass transport 
characteristics of polypyrrole films. For example, as mentioned above, the 
data of Buttry and coworkers indicates that Bu 4 N^ with a crystallographic 
radius tBu 4 N+° - 4.75 A (39) is excluded from pol 3 rpyrrole during 
electrochemical oxidation/reduction cycles whereas 0104 * ions (tcj^ 04 -° - 2.5 
A (39)) diffuse into and out of the film (30). In addition, scanning 
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electron micrographs (SEM) of polypyrrole obtained previously in this lab 
and by others (cf. 37,38) reveals that polypyrrole has an extremely dense, 
compact morphology. That individual pores or channels are not resolved by 
SEM indicates that any such structures must possess dimensions < ca. 50 A. 

However, several important discrepancies between the experimental and 
the simulated data obtained at 0.393 V (Fig. 4B) are evident. Most 
noticeable are differences in the transistion from the high frequency 45° 
data to the low frequency limiting region. In contrast to the simulated 
data, the experimentally observed transition is more gradual; occurring over 
a band of frequencies from 10 - 60 Hz. Keiser and coworkers have determined 
that the exact shape of the Nyquist plot for a porous electrode is sensitive 
to the geometry of the pores (25). Their calculations for a variety of pore 
geometries predict that positive deviations in the Nyquist data relative to 
the cylindrical pore case are observed in the transition region for pores 
which possess a decreasing radius (25); ie., a pore radius less at the 
electrode end than at the electrolyte end. A similar positive deviation is 
observed for the data obtained at 0.276 V (Fig. 4C) . This deviation may 
indicate that the mean pore diameter decreases as the polypyrrole film is 
traversed from the electrolyte/film interface to the electrode. 
Alternatively, this gradual transition may reflect the non-homogeneity of 
either r^ or 1. 

A second discrepancy relative to the simulated data is observed at 
frequencies above ca. 10 KHz, where Z[im] becomes negative. Experimental 
data in this frequency region was irreproducible from film to film. 

Negative Z[im] values at high frequencies reflect the presence of inductive 
effects which are not predicted by simple porous electrode theory (25,26). 
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Consequently, a more complex equivalent circuit may be required to 
thoroughly describe the electrical behavior of oxidized polypyrrole films. 

The Nyquist data obtained at 0.186 V (Fig. 4D) shows evidence of a high 
frequency semi-circle which is characteristic of charge transfer 
limitations. The emergence of this feature in the Nyquist plot is expected 
as the polymer begins to be reduced (19,29). The faradaic component of the 
current present at this potential may augment the purely nonfaradiac current 
due to the porous electrode character of the film and contribute to the 
shift in Ri which is observed (see intermediate data below). Thus, the 
equivalent circuit shown in Fig. 2B is not adequate to describe polypyrrole 
at this potential and the fit of the simulated to the experimental data in 
Fig. 4D is poor. 

Although no evidence for charge transfer limitations is observed at 
0.276 V (Fig. 4C) , the fit of this data to the cylindrical porous electrode 
model is intermediate to that of the 0.393 V data (Fig. 4B) and the 0.186 V 
data (Fig. 4D) . For example, transition to the vertical limiting region 
occurs at ca. 40 rad sec'^ (expt.) vs. 15 rad sec"^ (calc). At lower 
frequencies, the experimental Z[im] data is significantly higher than that 
predicted from Eqs . 15 and 16 . 

The fact that the fit between the experimental and simulated data 
improves as the potential increases from 0.186 V to 0.393 V suggests further 
improvement might be observed at yet higher open circuit potentials . 
Unfortunately, irreversible oxidation begins to occur at potentials greater 
than ca. +0.5 V vs. SCE and reliable AC data in this region cannot be 
obtained. However, the above data shows that the AC impedance of oxidized 
polypyrrole can be modeled reasonably well with an extremely simple finite, 
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porous electrode model which assumes a cylindrical pore geometry. Again, a 
better fit to the experimental data is possible if a decreasing radius pore 
geometry is considered. Nyquist diagrams for such films are easily 
calculated numerically using the procedure of Keiser et al. (25) . 

Lightly doped polypyrrole 

Figure 6 shows experimental and simulated Nyquist diagrams for 
polypyrrole at open circuit potentials of -0.333V, -0.35V, and -0.361 V. As 
the cyclic voltammogram in Figure 3 shows, these potentials are at the foot 
of the polymer oxidation wave. Since these potentials are well negative of 
£ 0 > (|EOf . Eq^I > 150 mV), the concentration of oxidized sites, Cppy + 

(where Cppy^ - Cgp 4 .), or the % doping of the polymer ought to be small. As 
discussed above, the data of previous workers (27,28) and the capacitance 
data of the present paper indicates that increases in the conductivity of 
polypyrrole at these Eqc values are small, ie. the polymer film is 
essentially nonconductive . Thus, the AC response of the system ought to 
conform to Limiting Case I. 

Figure 6A shows the experimental low frequency Nyquist plots at these 
three potentials. Note that Z[re] is nearly frequency independent for f < 5 
Kz and limiting behavior similar to that observed for oxidized polypyrrole 
is seen. Plots of Z[im] vs. at these frequencys are linear (r > 0.999) 

again showing that a constant, limiting low frequency capacitance obtains. 
Such limiting behavior is accounted for in Limiting Case I by a finite RC 
transmission line (Fig. 2A) where the nonfaradaic circuit components of the 
oxidized polymer circuit (Fig. 2B) are replaced by the charge transfer 
resistance, the Warburg resistance, R^, and the Warburg capacitance, 
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C^. In addition, the capacitance of the Pt substrate electrode, which was 
assumed to be negligible in the oxidized polymer case, must be accounted for 
in Limiting Case I since the total double -layer capacitance of the system is 
small . 

Figures 6B - 6D compare the experimental data obtained at higher 
frequencies with simulated Nyquist plots. Simulated data calculated from 
Eqs . 10 and 11 were fit to the experimental data using the diffusion 
coefficient, and the gradient of the potential with respect to the 

concentration, dE/dC 3 p 4 ^. The values of these parameters (see Table II) are 
identical to those calculated from the experimental R1 and Cl using Eqs. 13 
and 14. Note that the agreement between the experimental and the simulated 
data at all three potentials is good. The best- fit diffusion coefficients 
obtained from the simulation of ca. 1 x 10“^ cm^ sec“^ are approximately an 
order of magnitude faster than those measured by Genies and coworkers for 
polypyrrole (D - 2.3 x 10*^° cm^ sec'^ in Et 4 NBF 4 , MeCN) using 
chronoamperometry (40). As discussed above, this disparity is not 
unexpected in view of the very different nature of these two experiments; 
ie., in the potential step (Cottrell) experiment, diffusion coefficient 
information is obtained while the film is in flux from one redox state to 
another (29). However, the diffusion coefficients obtained here are in 
agreement with those obtained at identical films and potentials by small - 
amplitude current pulse experiments (41) . 

Diffusion coefficients can also be calculated from the semi- infinite 
<iiff'^sion region of the Nyquist data. Figure 7A shows the experimental 
Nyquist plots obtained in this region in greater detail. The slopes 
calculated from these data (Table II) are in close agreement with the 
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theoretically predicted slope of unity. Plots of the impedance modulus, 

|Z|, versus for this frequency domain are shown in Figure 7B. The 

diffusion coefficient can be calculated from these data by using values for 
dE/dCgp^. obtained from the finite diffusion region (Table II) and Eq. 12. 
The values for so obtained (Table II) are smaller at each potential 

than those calculated from the finite diffusion frequency domain; (Dsi/Dp). 
0.333 ” 0*28, (Dsi/Dp)- 0 . 3 5 0 ” 0.31, (Dsj/Dp).Q — .28. Note that Dsj/Dp 
is ca. 0.3 at all three potentials. Rubinstein et al. have observed similar 
behavior for the AC impedance of the reversible redox couple Ru(bpy) 3^'^^2 + 
in Nafion modified electrodes (20). At Nafion modified electrodes, however, 
diffusion coefficients obtained from the semi- infinite diffusion region are 
faster by a factor of ca. 3 than those obtained from the finite diffusion 
limiting region. 

The factors which might account for the observed discrepancy between 
Dgj and Dp in the polypyrrole films studied here are identical to those 
discussed in detail by Rubinstein and coworkers (20) . Note that an 
additional complication in the case of pol 3 rpyrrole is the estimation of 
dE/Ggp^. which cannot be determined by a straightforward coulometric 
titration without a reliable double-layer capacitance correction. 
Consequently, the value of dE/dCgp^. employed here is that derived from the 
finite diffusion data. The simultaneous determination of dE/dCgp^. and 
Dp, BP 4 . is susceptible to error from film thickness nonuniformity. However, 
because the thickness of the films used here is known to be 0.27 m® ± 0.02 
urn, film nonuniformity seems an unlikely source of the diffusion coefficient 
effect. 

The possibility exists, however, that the morphology of the polypyrrole 
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film is nonuniform over the film thickness. Under these circumstances, 
^Sl/^F < 1 would be produced by a film morphology which was dense 
(transport restrictive) at the electrode/film interface and relatively open 
at the film/electrolyte interface. It is interesting to note that the 
existence of such a morphological gradient could also explain the AC 
impedance data obtained for the oxidized polymer (see explanation above) . 
Furthermore, it is not unlikely that such gradients in the polymer void 
volume or effective pore diameter might be introduced during electrochemical 
synthesis of the polymer film. This is particularly true since the 
polypyrrole polymerization reaction is not diffusion controlled at the 
current densities employed here (1 mA cm'2)(40). Thus, pyrrole monomer is 
available at the film edge and the polymerization reaction can proceed from 
the interior surfaces of the film as well as at its exterior (apparent) 
surface. As a result, the total time available for polymerization at the 
interior surfaces of the film increases from the film/monomer interface to 
the electrode. The effective pore diameter of the electrochemically 
synthesized film, then, increases from the electrode surface to the film 
edge. It should be noted that although a morphological gradient in the 
polypyrrole film such as that described above is consistent with the 
observed AC impedance data, alternative explanations exist. For example, 
Rubinstein et al. have accounted for the similar effects observed at Naf ion- 
coated electrodes with a model which incorporates two uncoupled parallel 
diffusion pathways (20) . 

A second similarity between the present data and that of Rubinstein et 
al. is the behavior of the experimental data in the transition region (23). 
As shown in Figs. 4B - 4D, experimental Nyquist plots exhibit a gradual 
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transition from the semi- infinite (45°) to finite (vertical) diffusion 
regions of the Nyquist plot whereas a relatively abrupt transition is 
predicted by Eqs. 10 and 11. This effect is probably best attributed to 
inhomogeneities in the polymer structure which result in a distribution of 
diffusion pathways (apparent film thicknesses) . 

Intermediate oxidation levels. 

At potentials from -0.3 V to +0.1 V, significant concentrations of both 
oxidized and reduced polymer are present in a film which is also a good 
electronic conductor. Under these circumstances, neither of the equivalent 
circuits employed above are adequate to describe the AC impedance of the 
polymer. The faradaic impedance of porous electrodes has been discussed in 
detail by De Levie (32-34). The proposed models, however, do not consider 
the finite diffusion case and become complex if the effective pore diameter 
approaches dimensions of ca. 1 /jm (34). In addition, the modeling of 
polypyrrole at these intermediate potentials is complicated by the need to 
know the specific double layer capacitance (or the true surface area) of the 
polymer at each potential. More sophisticated porous electrode models (cf. 
42,43) introduce additional adjustable parameters which are difficult to 
experimentally identify. Finally, it should be noted that the work of 
Feldburg (35) and others suggests that significant extensions of established 
porous electrode theory are likely to be neccessary to fully account for the 
AC impedance behavior of conducting polymers. 

Although we have not attempted to identify an equivalent circuit 
appropriate for the AC data obtained at intermediate potentials, important 
qualitative data can be obtained in the absence of such a model. Figure 8 
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shows Nyquist plots for AC impedance data obtained at potentials from +0.1 V 
to -0.3 V vs. SCE. The low frequency data (Figs. 8A and 8B) show that these 
data are qualitatively similar to that which was observed at more positive 
and more negative potentials. Both the capacitance and the resistance 
assume limiting values, and , at frequencies less that 5 Hz. A plot of 
Rl and vs. Eqc obtained over the entire potential region from +0.4 V to - 
0.3 V is shown in Figure 9. Note in particular that the capacitance in the 
region from -0.3 V to -0.4 V is small relative to the maximum Cj_ values 
observed at higher potentials; ie., for the fully oxidized polymer. This 
data supports our assumption (Limiting Case I) that polypyrrole -BF4' is 
relatively nonelectronically conducting at these potentials. The value of 
Cl for the oxidized polymer (3 x 10-2 F cm'^ @ Eqc “ 0.393 V)) corresponds 
to a bulk capacitance of 350 F cm'^. This number is significantly higher 
than that observed by Burgemayer et al. (200 F cm'^ (2)) and by Diaz et al. 
(100 F cm" (36)). Note also that for a maximum is observed at ca. -0.2 
V. These data are similar to that obtained by Tanguy et al. for thick 
polypyrrole films (13,14). In the vicinity of the switching potential (E®') 
both the faradaic components (Ri^f and C]^ f) and the capacitive components 
(^l,c ^l,c) ^1 ^1 expected to be significant. 

To a first approximation, the limiting low frequency behavior of these 
components R]^ and can be inferred by superimposing the predictions of 
Limiting Case I (Eqs 13,14) and Limiting Case II (Eqs. 18,19). If the 
specific double layer capacitance, of polypyrrole can be assumed to 

possess an S-shaped dependence with potential (no maximum), Eq. 19 predicts 
that if r^ and 1 remain constant, vs. Eq^ will have a similar shape. 

In this case, since the faradaic component, C]^ f, is a maximum at Eq^ - 


25 



E°'(24), the superposition of faradaic and capacitive components yields a Cl 
vs. Eoc relation which has a maximum near Eo' , as observed here and 
previously (13,14). Equation 13 predicts that the faradaic component f 
ought to be a minimxam at Eq^ - E°'. Thus, it is difficult to account for 
the potential dependence of observed here. Since the capacitive 
component, R^ ,, is proportional to 1, and r^"^ (Eq. 19), relatively 

dramatic changes in the pore structure (eg. increases in rc with potential) 
are probably neccessary to account for the R^ vs. Eqc relationship obtained 
by us and by previous workers (13,14). 

Figure 8C,D shows the high frequency Nyquist data obtained at 
intermediate potentials. Of particular interest is the absence of a well 
defined semi circle resulting from charge transfer limitations at high 
frequencies (26,29). The AC impedence data obtained by various workers 
(2,13) for thick films invariably reveals the presence of such a semi-circle 
at high frequencies. Unfortunately, quantitative measurement of the 
heterogeneous rate constant for this system is impossible given the porous 
electrode nature of the film; ie. , the electroactive area of the film is not 
well defined. However, the present data suggests that the apparent rate 
constant may decrease considerably with increasing film thickness. 

The high frequency resistance, R^j, obtained at Z[im] — 0 is compared to 
C]^ in Figure 10. Rj^ values are good approximations of the total series 
resistance of the solution and the polypyrrole film (29) . The Rj^ values 
observed at positive potentials (R}^ - 6.5 ohms + 0.2 ohms) can be attributed 
entirely to the supporting electrolyte resistance; ie., the film resistance 
of the oxidized polymer is negligible. At potentials negative of -0.2 V, Rj^ 
increases rapidly to ca. 20 ohms at -0.6 V. It is interesting to note that 
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decreases in commensurate with oxidation of the polymer film occur in 
advance of increases in Ci. Note, for example, that at -0.3 V, the value of 
Rh is 15% of its value at -0.6 V whereas Cx has only attained ca. 15% of its 
limiting value. Since the value of Cx is sensitive only to the electronic 
conductivity of the film whereas that of R^ can reflect either ionic or 
electronic conductivity, these data suggest that polypyrrole is a good ionic 
conductor at potentials (near -0.3 V) where it is a relatively poor 
electronic conductor. 

Conclus ions 

AC impedance data has been presented for thin polypyrrole film- 
modified electrode surfaces . Although an equivalent circuit which accounts 
for the observed AC behavior of polypyrrole over a wide range of potentials 
is likely to be complex, two limiting cases exist for which the equivalent 
circuits are simplified. At low doping levels (Limiting Case I), the 
equivalent circuit appropriate for polypyrrole is identical to that used 
previously to describe the AC impedance of Nafion polymer modified 
electrodes (20.24) and metal oxide, Li intercalation films (19,21-23). 
Diffusion coefficients were obtained from the analysis of Ho et al. (19) for 
both semi- infinite and finite time domains. The approximate ratio of the 
values Dgx/Dp ” 0.3 at the potentials examined. This observation is 
consistant with the existence of a morphology gradient normal to the surface 
of the electrode where the effective pore diameter of the pol)nner, and 
therefore the diffusion coefficient Dgp^. , increases in the film with 
distance from the electrode surface . 

The AC impedance data observed for the oxidized pol 3 mer (Limiting Case 
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II) was compared to that predicted for finite, porous metallic films 
(25,26). The agreement between the experimental and the simulated Nyquist 
data (calculated from ref. 26) improves with increasing Eqc- Using the 
analysis of Candy and coworkers, film parameters such as the pore radius, 
pore length, and pore number density can be expressed as their cylindrical 
pore electrode equivalents (26). The values so obtained for thin 
polypyrrole films are consistent with what is currently known about mass 
transport in such films. Better agreement between the simulated and the 
experimental Nyquist plots ought to be obtained if other pore geometries are 
considered. In particular, the transition region of the Nyquist data is 
better accomodated by a model in which the pores possess a decreasing radius 
from the film edge to the electrode surface. This conclusion is in 
qualitative agreement with the diffusion coefficient data obtained for the 
oxidized polymer. 

Thus, both mass transport and structural information can be gleaned 
from AC impedance analysis of polypyrrole films based on these two 
relatively simple equivalent circuits. In addition, the potential 
dependence of the low frequency series resistance and capacitance were 
determined and the origins of the observed dependencies were discussed. 
Comparisons between the high frequency resistance and the low frequency 
capacitance Indicate that significant ionic conductivity is present for 
polypyrrole films at very low oxidation levels, prior to the onset of 
significant electronic conductivity. 
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Table I. Sunnaary of Film Parameters Calculated From 
AC Impedance Data For Oxidized Polypyrrole . 


equivalent cylindrical pore number 


Eoc 

pore radius® 

density^ 

pore length 

. ys, sgg) 

(cm) 


( cm) 

0.393 

4.26 X 10*8 

5.43 X 10^2 

8.72 X 10*^ 

0.276 

2.63 X 10*8 

1.23 X 10l3 

1.01 X 10*^ 

0.186 

1.96 X 10*8 

1.90 X 10^3 

1.17 X 10*2 


^Values for r^, and n were determined by fitting Nyquist plots calculated from 
Eqs. 15 and 16 to the experiment data. The values of C]^ used for these 
calculations were obtained from the slope of Z[im] - (wCg)'^ vs. at f < 5 Hz. 

was determined from Z[re] at ca. f - 1 Hz. and corrected for R^. 

^Calculated from Eq. 20 assuming the pore volume Vp - the total film volume, 



Table II. Summary of Transport Data Calculated For Semi-Infinite (SI) and Finite (F) 

Frequency Regimes For Lightly Doped Polypyrrole. 


Eqc 

Dp^ 

dE/dC^ 

d|Z|/d(w‘^/2)b 

DSI^ 

(V vs. SCE') 

(cm^ sec-=-^) 

(V cm^ mol-^') 

(ohms rad-^ sec) 

(cm2 sec-1) 

-0.333 

1.40 X 10’^ 

593 

1.36 

o 

1 

o 
»— 1 

00 

-0.350 

3.31 X 10-9 

1140 

1.02 

1.03 X 10*9 

-0.361 

3.38 X 10-9 

1300 

1.20 

9.36 X 10-^0 


^Obtained by fitting Nyquist data calculated from Eqs. 10 and 11 to the experiment 
Nyquist data at each potential. The best fit values Dp and dE/dC so obtained approximate 
those obtained from R]^ and and Eqs. 13 and 14. 

^Obtained from the slope of the 45° linear region of the Nyquist plots shown in 
Fig. 7A. 

^Calculated from d| Z| /d(w* ^ ' 2 ) using Eq. 12. 



Figure Captions 


Figure 1 - 4 ^ Randel's equivalent circuit with a finite Warburg element, 

(^9) > Sj. Finite RC transmission line circuit equivalent to 
that of a finite pore with a uniform radius (ex. a cylindrical 
pore) . 

Figure 2 - Schematic diagram of the platinum disk electrode employed in 

this work. a. 7mm glass tube, b. copper wire, c. Kel-F electrode 
body; diameter - ca. 10 nun, d. silver epoxy contact, e. Pt disk; 
diameter - 3.3 mm. 

Figure 3 - Cyclic voltammogram for a typical 0.27 fim polypyrrole film in 
0.2 M BU 4 NBF 4 , MeCN electrolyte at 20 mV sec'^. 

Figure 4 - 4 ^ Nyquist diagram comparing the low frequency (f < 0.25 Hz.) 

AC impedance data obtained at three potentials for an oxidized 
polypyrrole film. C - 0.393 V, • - 0.286 V, A - 0.186 V. 
Frequencies in Hz. are as shown. 

Bi£j. Nyquist diagrams comparing experimental and simulated 
AC impedance data for an oxidized polypyrrole film. Frequencies 
(Hz.) for simulated data (solid line) are indexed on the right of 
each plot. Frequencies for the experimental data (open circles) 
are indexed at the left of each plot. 

A^ 0.393 V, ^ 0.286 V, 0.186 v! 

Simulated data calculated from Eqs . 15 and 16. 



Figure 5 - Typical plot of Z[im] vs. u>'^ for data obtained at low 

frequencies. These data at f < 5 Hz. and 0.393 V. 

Figure 6 - ^ Nyquist diagrams comparing the low frequency A.C impedance 

behavior at three potentials for a lightly doped polypyrrole 
film. O - -0.333 V, A - -0.350 V, • - -0.361 V. 

B-D . Comparisons of experimental and simulated Nyquist diagrams. 
Frequencies (Hz.) for simulated data (solid line) are indexed on 
the right of each plot. Experimental data (open circles) are 
indexed at the left of each plot. -0.333 V, ^ -0.350 .V, 

-0.361 V. 

Figure 7 - ^ Comparison of the 45° semi- infinite region of the Nyquist data 
obtained for a lightly doped polypyrrole film. The slopes 

calculated for these plots are listed in Table II . 

O - -0.333 V, • - -0.350 V. A - -0.361 V. 

B. Plots of |Z| vs. for the linear regions of Che plots 

shown in Fig. 7A. 

□ - -0.333 V, A - -0.350 V, O - -0.361 V. 

Figure 8 - A.B. Comparison of experimental Nyquist data obtained at low 

frequencies for potentials from +0.1 V to -0.3 V. Frequencies 
in Hz. as shown. 

4^ 0 - 0.093 V, • - 0.001 V, A - -0.091 V, 

A - -0.140 V. 

^ O - -0.211 V, • - -0.255 V A - -0.299 V. 



Figure 8 - C.D. Comparison of higher frequency experimental Nyquist data 

obtained at the same potentials. Identity of plots same as for 
Fig. 8 A,B. Frequencies in Hz. as shown. 

Figure 9 - Low frequency resistance, and capacitance, C]^, vs. Eqc 
measured at 0.25 Hz., obtained from plots of Z[im] vs. 

Figure 10 - High frequency resistance, R^, and low frequency capacitance. 
Cl, vs. Eoc. R^ - Z[re] at Z(im] - 0. 
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